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For a one-dimensional electron-phonon system we consider the photon absorption involving elec- 
tronic excitations within the pseudogap energy range. Within the adiabatic approximation for the 
electron - phonon interactions these processes are described by nonlinear configurations of an instan- 
ton type. We calculate intensities of the photoelectron spectroscopy PES including the momentum 
resolved one ARPES and supplement to known results for the optical subgap absorption. We start 
with the generic case of a one dimensional semiconductor with pronounced polaronic effect. In 
details we consider the Peierls model for a half-filled band of electrons coupled to the lattice which 
describes the polyacethylene and some commensurate Charge Density Waves. Particular attention 
was required for studies of momentum dependencies for the ARPES where we face an intriguing 
interference between the time evolution and the translational motion of the instantons. 

PACS numbers: 72.15.Nj 78.40. Me 78.70. Dm 71.45.Lr 



I. INTRODUCTION. 

This article is devoted to theory of pseudogaps in electronic spectra as they can be observed by means of the Photo 
Electron Spectrography (PES) or the Angle Resolved Photo Electron Spectroscopy (ARPES). Notion of a pseudogap 
(PG) refers to various systems where a gap E g in their bare electronic spectra is partly-filled showing subgap tails. 
The best known examples are the tails in the Density jpf States (DOS) due to disorder or the Urbach tails in the 
subgap optical absorption due to thermal fluctuations! 2 !. But the PG is especially pronounced in cases where the 
bare gWWS opened spontaneously as a symmetry breaking effect which was the subject of detailed experimental 
studiesErQ. In quasi-l£> conductors it is known as the Peierls-Frohlich instability leading to the Charge Density Wave 
(CDW) formation, as well as to-analogous spin-Peierls and Spin Density Wave statesB. Here the picture of the PG 
has been suggested theoretically^ in relation to absence of the long range order (LRO) in ID CDWs. Similarly in 
High-Tc materials the gap opening mechanisms (still of a disputable origin) which are not stabilized yet by the LRO 
are considered. They are supposed to be responsible for the pseudogap which opens_at higher energy scale than 
the sharp gap appearing below the transition temperature to the superconducting stateliJ. We shall consider generic 
ID semiconductors and concentrate on systems with the dimerized ground state like the well known polyacethylene 
(Cif ^Mt 2 ] or some CDWs with the 2-fold commensurabilitv like NbSjB. Properties of incommensurate CDWs (blue 
bronzes, tri- and tetra-chalcogenides of transition metals )Ij arc further complicated by interference of the gapless 
collective mode and we shall consider them separately. 

For ID systems with only a discrete symmetry an absence of the LRO is not drastically important at low enough 
temperatures. At finite temperature there is a remnant concentration n s ~ exp[— E s /T] of topological solitons (kinks 
with the energy E s ) commuting different domains of the order parameter ~ ±Ao. The midgap electronic states 
associated witlu-these solitons create a finite DOS ~ n s at the former Fermi level E = and originate related optical 
features (seeEj'E£l for a review andlij for a systematic theory), hence there is no true gap at finite T. Neglecting this 
exponentially freezing contribution, we are left, at first sight, with a sharp electronic gap E® = Aq. But what happens 
instead is that, even at T = 0, rather smeared edges appear at ±Ao while the spectrum extends deeply inwards the 
gap. This effect is particularly pronounced in ID because of the edge singularity ~ (E — Eg) -1 ' 2 at E as E® in 
the bare DOS which smearing is considered to be the most evident signature of the PG formation below E g . In the 
Peierls state the PG effect is further enhanced by strong electron-phonon (e-ph) interactions between electrons and 
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quantum fluctuations of the gap amplitude 5 = A(x,t) — Aq. Stationary excitations (eigenstates of the total e-ph 
systemXare now the selftrapped states, polarons or solitons, which energies W p or W s are fractions of A , see the 
reviewed. The states close to the bare electron edge Ao can be observed only via instantaneous measurements like 
optical or X-ray absorptions or tunneling. The spectra of these nonstationary states fill the range Ao > E > W p for 
single electronic (PES) or 2Ao > u> > 2W S for e — h (optics) processes, particularly near Ao the states resemble free 
electrons in the field of uncorrelated quantum fluctuations of the latticeEj. Here the selftrapping has not enough time 
to be developed. But approaching the exact threshold they evolve towards eigenstate which are selftrapped states 
accompanied by excitations of their dress. This picture describes coexistence of the PG region Ao > \E\ > W P ,W S 
and the exact gap \E\ < W p , W s in a similarity with the High-T c superconductors. 

It should be stressed in this respect that there cannot be a common PG for processes .characterized by different time 
scales. We should distinguish ( as stated ir£3 and stressed also in several reviews, e.g.Eil) between short living states 
observed in optical, EPS (and may be tunneling) experiments and long living states (amplitude solitons, phase solitons) 
contributing to the spin susceptibility, NMR relaxation, heat capacitance, conductivity, etc. States forming the optical 
PG are created instantaneously: over times which are shorter than the inverse phonon frequencies r opt ~ n/E g < <jSZ h 
and many orders of magnitude beyond the life times required for current carriers, and even much longer times for 
thermodynamic contributions. It follows then that the analysis of different groups of experimental data within the 
same modelO must be reevaluated. 

While experimental techniques of PES and ARPES are still less accurate than the traditional optics, their very 
fast progress through the last decade allows to rely upon a required accuracy already today or in nearest future. 
An apparent advantage of the ARPES is its access to the momentum distributions of spectral densities. Another 
particular feature of both PES and its ARPES version in our applications is its decoupling from the final state 
Coulomb interactions which affect drastically the optical intergap absorption. In this respect the PES differs also 
from the traditional spectroscopy by internal absorption from a core level; their the final state interaction batween 
the band particle and the remnant hole, originates a seminal problem of the X-ray edge singularity in metalsa. 

A theory of-the subgap absorption in optics has been developed already: for a general type of polaronic 
semiconductorsH with an emphasis to long-, range Coulomb effects and for the one dimensional Peierls type system 
with an emphasis to solitonic processesOl23 - the last work is closest to our targets. Here we shall address the PES 
and ARPES in ID systems. Methodologically all these studies deal with quantum transitions between distant field 
configurations which generalize the WKB method for a single degree of freedom. These classically forbidden nonlinear 
processes, well localized in space and time, are described by extremal trajectories called the instantons. The methods 
of op±ixnal trajectories in the functipaal space have been initiated by the problem of quantum decay of a metastable 
stateEa (that is of |*u" false vacuum" EJ) and for quantization of solitonS (see the bookEj for a review). The relevant 
article on this linecil was devoted to selftrapping barriers for transformation of a bare electron to the polaron in 3D 
systems. 

The plan of our article is the following. In Sec. II we describe a general scheme for calculating the ARPES 
probabilities within the adiabatic approximations, that is based on smallness of typical phonon frequencies with 
respect to energies of electron transitions. In Sec. Ill we explore applications to particular models. In Sec. IIIA 
we consider a tutorial zero- dimensional model for an electron coupled to a single oscillator mode and show that 
results of the adiabatic approximation coincide with exact calculations. In Sec. IIIB we study the transition rate 
for shallow subgap states in a ID system, both near the bare gap and at the polaronic threshold. In Sec. IIIC we 
consider in details the Peierls model in the half-filled band case. In Sec. ItfA we present the ARPES theory for 
this model. Beyond the PES, in Sec. HIE we also extend earlier calculationsc3E9 for the intragap optical absorption 
related to creation of kinks pairs. In Sec. IIIF we consider relations of the general approach to the model of quantum 
instantaneous disorder. Finally Sec. IV is devoted to discussions conclusions. 



II. GENERAL RELATIONS. 



A. Adiabatic approximation. 

The ARPEsi means an absorption of a high energy Slo photon taking off an electron from the crystal which is then 
analyzed in energy E out and momentum P out . Thus one obtains an information on the spectral density of the hole 
left beyond. The transition rate 1(^2, P), where Q = £1q — E out and P = —P ou t, is proportional to the imaginary part 
of the single-electron retarded Green function 

I(P,Q) oc/m J dxe- iPx J dTe ira G(:r,T,0,0). (1) 
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(Since now on we shall omit all constant factors and take the Plank constant % = 1; fi will be measured with respect 
to a convenient level: the band edge or the middle of the gap.) The simple PES nonresolved in momenta measures 
the integrated absorption intensity = J dpl(p, f2). 

We shall use the adiabatic (Born-Oppenheimer) approximation. Electrons are moving in the slowly varying phonon 
potential, e.g. Q(x,t), so that at any instance t their energies Ej(t) and wave functions ipj(x,t) are defined from 
the stationary Schroedinger equation for the instantaneous lattice configuration and they depend on time only para- 
metrically. The intensity can be written in the form of a functional integral D[Q(x,t)] over lattice configurations 
Q(x,t) 

/(0,P)oc J dxe- iPx J dT J D[Q(x,t)}^ o (x,T;[Q})%(0,Q;[Q})exp(-S[Q}) (2) 

This equat ion is already written in the Euclidean space it —> t which is adequate for studies of classically forbidden 
processesBESJ. Here ^o(x,t; [Q]) is the wave function of the added (N + l)th particle in the instantaneous field 
Q(x,t) which corresponds to the energy level E (t) = E [Q(x,t)] inside the gap. The effective action S = S[Q(x,t)] 
is expressed via Lagrangians Lj as 

M r-T r-oa 

S= dtL + / dt(L t - Q) + / dtL (3) 

Joo JO JT 

where indices j = 1,0 label systems with and without the additional particle state (actually the hole). Their typical 
structure is 



Lj = j dx ~{d t QY + Vj(Q) , n = const ~ w^ 2 (4) 

where uiq is a bare phonon frequency (implying the dispersionless phonon branch). A nature of the field Q, its 
potentials Vj(Q) and the action upon electrons differ for various models as described in Sec. III. Usually 
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Vi [Q(x, t)} = J dx^ + ej [Q] , ag 2 = uj^ (5) 

which contains the bare harmonic term ~ Q 2 and the adiabatic contribution from the energy of the electron system 
in the j-th excited state. The electronic terms j = 0, 1 correspond to adiabatic ground states of N and N — lelectrons 
in the instantaneous field Q(x). For calculations of subgap processes only lowest localized states j = 0, 1 are relevant, 
while other states belong to the continuum spectrum above the gap. Then ^ in (||) is the wave function ^ of the 
split off singly filled electronic state with the energy level Eq inside the gap and V\ — Vq + Eq, 

Within the PG region the main contribution to I comes from the saddle points of the action S: SS/SQ = 
dS/dT — 0, the last equation determines the value T = T(f2) as 

L^T) - L (T) = E (T) = Q. (6) 

(The same relation holds at the point t = if one substitute 0, T — > t, t + T in Eq. (0) and find the minimum over 
t). We shall find also, in Sec. HID devoted to particularities of the ARPES, circumstances when the extremum must 
be determined for the whole expression under the integral in the definition (Eh, including the wave functions in the 
prefactor of (^|). Any nonzero contribution requires for a finite action So < oo which selects configurations Q(x,t) 
deviating from the ground state only within a finite space-time region. Such extremal solutions with finite actions are 
called the instantonscJ, their trajectories correspond to tunneling in the real time. 



B. Translational mode. 



Within the saddle point approximation the transition rate is given by 

7(p,Q)=7 exp(-5o), (7) 

The prefactor Iq = ItIx ■ ■ ■ comes from integration over deviations SQ, 8T around the extremal solution. E.g. the 

factor I T = (d 2 S/0T 2 y 1/2 = y/dT/dCl comes from the integration over T. Usually it can be taken after the extremal 
is determined. But there are cases when J should be determined selfconsistently, like in statistical physics looking for 
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the minimum of a free energy, rather then simply of an energy. This option appears naturally in comparison of PES 
and ARPES where it comes from treatment of the translational invariance. Apparently an essential contribution to the 
action from the "quantum entropy" SS = — In Jo may come only from integration over particular zero modes appearing 
because of continuous degeneracy. The contribution of usual nondegenerate modes must be small by definition of a 
well defined extremum. 

The extremal solution can be written as Qo(x — X, t), x So(x — X, t) where the collective coordinate X of the instanton 
corresponds to the translational invariance which originates one zero mode in the functional integration around the 
saddle point solution. We expand the field Q(x, t) in the vicinity of instanton solution as 



Q(x, t) = Q (x - X(t), t) + r]{x ~ X(t), t) 



(8) 



with the orthogonality condition J dxQoi] = 0, wJaere r/(x, t) contains only nonzero modes. Integrating over the zero 
mode is carried out by means of Faddeev-PopovEj method inserting the identity 



1 = J D[X(t)]5(J dxr) {x-X(t),t)Q(x,t))J, J = J| / dx9 *Vo{x - X(t),t)Q(x,t), 



(9) 



where 770 = d x Qo(x — X(t) 7 t)/y J(d x Q) 2 dx is the normalized zero mode eigenfunction. The integration over X(t) 
can be done exactly. Taking into account only terms containing X(t) we rewrite (0) as 



I 



x oc 



dxe 



D[X(t)]J* (x - X(T),T; [Q])%(-X(0), 0; [Q]) exp 



~ / dtM(t)X 2 (t) 



where J = Ylt \J M(ti), and M(t) is the translational effective mass of the instanton at a given moment t: 

M(t) = n [ dx(d x Q) 2 . 



Introducing the integration over X\ = X(0), X2 = X(T) one obtains 
dXi dX 2 



ix oc 



dxe 



-iPx 



D[Q(x,t)]V (x - X(T),T; [Q])%(-X (0), 0; [Q])hhh, 



where 



r X(Q)=X 1 rO ,,,, >2 , r X(T)=X 2 r-T 

h= Dm)] n ^MW)e- f -~ dtMX/ \h= Dm)} n vme-^ dtMx2 '\ 

JX(-oo)=0 „„, JX(0)=X! ue(o,T) 



ti£(-oo,0) 



(10) 



(11) 



h= [ XiOO)=0 DlXit)] [] ^MttT)e-fr dtMA2/2 . 

Jx(T)=x 2 



Each integral li is easy calculated, after the transformation MX 2 = Y 2 , with the help of the resul 



E 



x(T)=x 2 



nI /,m - [ T (x 2 + V(x))dt 



d 2 S 



dx\dx2 



1/2 



-So 



' x(0)— x\ 

where is the saddle point action. After simple integrations over dx, dXi, dX 2 we arrive at 

fT 



Ix oc| * P (T) \ 2 exp 



P 2 

„ 2Mit) 



(12) 



where ^ p is the Fourier transform of the wave function 

Contrary to stationary solutions, the solitons (polarons, kinks), here M(t) depends on time along the instanton 
trajectory Qoix, t). The exponential term in (12) gives the contribution to the total action which is important at very 
large P because generally it is of the order of M -1 ~ cOq. It provides an important, sometimes leading effect upon 
the edge shape: it is singular along the instanton term with M vanishing at large \t\ when the space localization is 
week and Af _1 (t) diverges. 
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C. Zero dimensional reduction. 



As a rule the space-time differential nonlinear equations for extremals are not solvable, hence one is bound to 
variational procedures. We shall use the following ansatz which is actually an one-parameter reduction of the functional 
space. The instanton trajectory must satisfy the condition (^) which means that there is an electron level Eo = Q 
at times t = 0, T. It seems reasonable that the trajectory with only one local level inside the gap will be a good 
approximation for the extremal action. There is also a special advantage that for models considered below we know 
exact solutions of the stationary problem Qo(x,a) depending on some continuous parameter a with the local level 
E (a) inside the gap. Therefore we shall search the instanton (time-dependent) solution in the form Qo(x, a(t)). After 
integration around saddle point Da(t) we find finally for ARPES intensity 

I(P,n)^^\^ p (a T )\ 2 e- s °, (13) 



where So is the instanton action which is extremum over T, a(t) of 

P 2 
2M{a{t)) 

where 



roc pi 
S= dt[f(a)d 2 + ——-- + V(a(t)}, (14) 



ft \ V I A ( dQa(x,t)\ 2 

fia) = 2j dx \-d^T) (15) 

is the variable effective mass, V(a(t)) = Vi(a(t)) — ft, for £ [0, T], or V(a(t)) = Vo(a(t)) otherwise. 

In some cases the form- factor "fp gives the main contribution to the function /(P, ft). Then we must find the 
extremal solution for a modified action S — > S — log | "fp | 2 . Such a possibility will be considered in details in Sec. 
HID. 

Since the PES intensity is obtained simply by integration of ARPES one, we obtain 

' '' '111, > 

where 




m <x I / g e -P < | Vp(T) I" ] \j — exp[-S ], (16) 



I 2 = 



dt 



o M(a(t)Y 



and So is the extremum of the action ( |l4| ) without the kinetic term P 2 /2M. The form-factor F = J dP \ tI> p (T) | 2 
exp(— P 2 l 2 /2) in (|l^ ) provides a kind of Debye- Waller reduction of intensity. Thus for small / in compare to the 
localization length £ of ^ (short times near the free edge) F m 1 while for large I ^ £ (long times T near the absolute 
threshold) F ~ \^ p = q(T)\ 2 ~ £/Z with £ being the polaronic width. In our examples typically it will hold I <C £. 

For typical potentials V(a) (see Fig. 1-4 below) the extremum solution with a finite action exists only in some 
region of f2 where the potential curve Vi(a) — ft (or Vi(a) — ft + P 2 /2M(a) for the case of ARPES) crosses both the 
curve Vo(a) and V = 0. It takes place if the minimum of the potential V\ — O ( or Vi — O + P 2 /2M) is placed below 
the minimum of the potential Vq (min Vq — 0). If ap is the point where 

Vb(a T ) = Vi(a T ) - ft (PES), V (a T ) = Vi(a T ) - ft + P 2 /2M{a T ) (ARPES), (17) 
and a/ is the solution of 

T4(a/)-ft = (PES), Vi( a/ ) - ft + P 2 / 2M ( a /) = (ARPES), (18) 
then the extremal solution is the closed trajectory, described by the equation 

f(a)d 2 = V-(o), (19) 

where V — V (a), a € [0, a T ], and V — Vi(a) — ft + P 2 /2M(a)}, a e (aT,a/] for ARPES, or the same expression 
without the kinetic potential P 2 /2M for PES. The Eq. @ is a consequence of the Eq. (|6|), it implies a continuity 
of the velocity a at the point qt- The Eq. ( fl8| ) reads that the velocity d = at the point aj. 
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The instanton equation (|l9|) describes the motion of some particle of zero energy with variable mass 2/ (a) in the 
inverted potential —V. The trajectory starts at t = — oo from the point a = 0, reaches the points ax at t — 0, a/ at 
t = T/2, after which the particle moves back across the point ax at t = T to the initial point a = at t — > oo. The 
instanton action (for ARPES) can be expressed as 

fir rqf 

So = 4/ dq^/fVo + 4 <W/(^i-^ + -P72A/(a)). (20) 

J J qt 



III. RESULTS FOR ABSORPTION INTENSITIES. 



A. Zero-dimensional case 



We shall start with a tutorial example where the general approximate scheme can be compared with exact calcu- 
lations. Consider a particle interacting with a single quantum oscillator at zero temperature. Physically it can be a 
problem of the Jahn- Teller center or a zero order approximation for the small radius polaron. The total Lagrangian 
is 



L 



-T 



^lq 2 - gqN, 



(21) 



where q is the oscillator coordinate, m is its mass, g is the interaction constant, and N is the electron's occupation 
number. The potentials Vq and V\ 



V 



rauj, 



V, 



gq 



muj. 



^(q-q ) 2 -W ; q 



g 



Wo 
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are shown in Fig. 1. In compare to TV = 0, for N — 1 the oscillator states n are shifted by qo in coordinate and by 
—Wo in energy. Adiabatically allowed transitions take place only at the given q that is they explore the region q = 0, 
il = 0. Taking into account the quantum character of the coordinate q, the transitions are found to the lower energies 
down to Wo- The region > fl > —Wo corresponds to pseudogaps in more complex systems. 




FIG. 1. Potentials V , Vi - O. 



The solutions of (|17|), (gj) read 

q T = -Q,/g, q f = 2g/(mw%). 
The instanton action is easy calculated and for the absorption rate we arrive at 



G 



7(0) cx 



1 

y/il + W 



exp( ) exp 



"O + Wp 



log 



eW 

(o + Wo). 



Expanding the exponent near the extremum point £1 = 

i~(0) oc exp 



we obtain 

n 2 



2W cu 



(22) 



(23) 



which describes smearing of the adiabatically allowed electronic edge. This is the Gaussian with the width SSI = 
i/WolJo 3> which justifies the adiabatic approximation. At the lower boundary SI = —Wq of the exact spectrum 
the probability has a finite exponentially small value which is approached with an infinite slope. 

For this simple example W(Sl) can be calculated exactly in the representations of eigenstates. Indeed, 



I (SI) = ^|<n, 1|*+|0,0)| 2 ^ - E n + E ), 



(24) 



where |0,0) is the ground-state of the system without electron with the energy E — Hujq/2, and |n, 1) is the n-th 
excited state (n = 0, 1, ...) with one electron and the energy E n = hoj (n + 1/2) — W - The wave functions are the 
ones for the harmonic oscillator which are centered at q = for |0,0) and at g = go for \n, 1). In the limit ujq — ► 
and n = (0 + Wq)/u>o > 1 we find 



1(0) = 



exp[ 



-(<52/) 2 /2l 



■ exp 



nlog 



e{8y) 2 
2n 



(25) 



where Sy = g/ \J mwjj. This equation reproduces the result ( p2| 



B. Shallow states in a ID system. 



Consider the electron (hole) states in a ID dielectric near an edge of a conducting (valence) band. We shall take into 
account the dispersionless phonon field Q(x,t) with the bare frequency loq which interacts locally with the electron 
via the deformation potential with a coupling g. Within the adiabatic approximation the action has the form (again 
at imaginary time) 



S = 




dx 



It is well known that the stationary (Q = 0) extremum of S corresponds to the selftrapped state (the polaron) 
arrive directly at Eqs. ||,H with E and ^ as solutions of the Schrodinger equation: 



1 d 2 ^ 



gQ^> = m, E = E(t) = E[Q(x, t] t 



(26) 



We 



(27) 



2m dx 2 

The minimum of the initial potential min Vq = is achieved at the uniform configuration Q = 0. The minimum of V\ 
takes place at the polaron state minVi = W p = — Wo — —mg 2 /2A. The minimum of V\ — O must be lower than the 
minimum of Vq, hence —Wo < O < 0. 

Since the exact solution of the complete time dependent extremal equations is not known, we shall follow a variational 
approach. We will search for the instanton solution in the form which simulates the stationary solution but with the 
variable time dependent parameter B(t) which we chose as the inverse localization length: 



Q(x,t) = - 



B 2 (t) 



1 



mg cosh B(t)x 



* = 



B{t) 



1 



2 coshB^x 1 



E(t) = 



2m 



(28) 



The stationary solution corresponds to B(t) = Bq = mg 2 /2. In terms of b(t) — B(t)/Bo the Lagrangians become 



in = Wn 



%b 2 



2¥ 



%b 2 



2b A - W 



(29) 



where Co = 4 + 2-7r 2 /15 w 5.3. The instanton equation describes the motion of a zero energy particle, with a coordinate 
6 and with a variable mass / ~ 6, in the inverted potentials 

V = W 2b 3 and V x - O = W [2b 3 - 3b 2 ] 

which are shown in Fig. 2. 



O 
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The value &t = vKV^T' the vame bf at the turning point is determined by the condition Vi(bf) — ft = 0. 
For the PES the action is 



Wo 
u 



V2b 2 db 



{2b 3 -3b 2 -^-)bdb 



and we arrive at the following results. 

Near the free electronic edge ft = we have 

/(ft) oc| ft r 1 / 4 



exp 



^o~ (-ft) 3/2 
6 ojoWI 12 



(30) 



The characteristic width of the edge is (5ft ~ (woW-o) 1 / 3 3> Wo which justifies the adiabatic approximation. 
In the vicinity of the absolute edge ft « —Wo we obtain 



i (ft cx = exp —2 — exp 

V ' V^ + Wo 1 loo 1 



} C n + W ^e(V3-l) 2 Wo 



to 



(ft + Wo 



(31) 



In this limit the probability is finite, exponentially small in the adiabatic parameter Wo /too and h shows a weak 
singularity in its derivative. The prefactors in above formulas come from It ~ (d 2 S/dfl 2 ) 1 ' 2 . The form-factor is 



F 1 according to Sec. II. Indeed, in spite of a formal divergence of I, within limits of the adiabatic approximation 
(|ft| ^> too or (ft + Wo) ^S> too) it stays below the localization length £. Thus l/£ ~ (wo/ft) 1 ^ 2 *C 1 for small ft and 
l/t ~ (to m[W /(ft + Wo)]) 1 ' 2 for ft « -W . 

Peculiarities of the ARPES will be studied in the next section within a richer Peierls model. 



C. The Peierls model 



We consider now the PES absorption spectrum in the gap region for a half-filled Peierls model. In the ground 
state the electron spectrum has the form E 2 = Vpp 2 + Aq and since now on we shall put the Fermi velocity Vp = 1. 
The gap 2A ~ e^expf— 1/A] is opened as a result of symmetry breaking lattice deformations. We shall consider 
the case of a dimerization (trans-Polyacethylene, Spin-Peierls systems) which is described by the model with the real 
order parameter A taking equilibrium values ±Ao- The excited states are solitons (kinks), polarons and bisplitons 
(kink-antikink pairs) which are characterized by electron levels localized deeply within the gap (see the reviewEfl) . The 
adiabatic approximation is valid when the electron transition energies are much larger than the phonon frequency 
loo *C \Ei — Ej\. For characteristic \Ei — Ej\ ~ Ao this constraint coincides with the applicability condition for 
the Peierls model in generate. Away from the ground state the field A(x, t) extends a role of the field Q from the 
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previous section. The first difference is that A^O already in the ground state. The second one is that the whole sea 
of electrons at E < — Aq contributes to the selftrapping energyE.3 and to the instanton terms. 

The effective Euclidean action S consists of the kinetic and the bare potential lattice energy and of the sum over 
filled electron levels E a . 



S{A(x,t)} = [ dxLjdt, Lj = ( dx f » + Vj[A(x,t)} 
J J ttAu;^ 



where u> p h — VA^o is the bare frequency of the phonon responsible for the dimerization while luq is the amplitude 
mode frequency in the dimerized state A / 0. The index j characterizes perturbations in the set {<x,} of filled 
electronic states in the field A. The potential energy in the j-th electronic state is 



V$[A(a:,t)] = J dx^+Y,E ai {A(x,t)} -W gs , 



which we have defined relative to the ground state energy W gs at A(x,t) = Ao for a system of N = N a electrons (N a 
is the number of sites in the chain). 

The PES absorption adds one particle (hole) to the system, then the polaron state is formed by local deformations 
of the field— A (a;) which originate the pair of split off electronic levels ±Eq localized deeply in the gap. In optical 
absorptionEj an electron is excited across the gap, then a pair of distant kinks is formed with electron level placed 
exactly in the center of the gapJ-Ep = 0). The crossover solution describing the states with just one pair of localized 
electronic levels ±Eq is knownESEj: 

A s (x) = Aq(1 — tanha[tanh(A a;tanha + — ) — tanh(A :ctanha — ^)])- (32) 

Depending on the parameter a it describes evolution from the shallow polaron at a — > to the pair of kinks at a — ► oo. 
In the first case the parameter a -C 1 becomes equivalent to b from the previous section. In the second case a I 
becomes a distance between divergent kinks. For the configuration |3^ the potential energy functional V and the 
energy of the local level Eq are given by 

E _ A 



K(a) = vE + - V A o - E o - ~ E o cos" 1 , E = (33) 
7r v 7r Ao cosh a 

where v = 0, 1,2 is the number of particles (electrons and holes) added to the system. The wave function of the 
intragap state can be written as0 ^o(^) y/ ~ A s (x) 2 . 

The case v = describes the state which evolves from the ground one without perturbing the occupation numbers: 
the level — Eq is doubly filled while the level Eq is empty. Naturally the minimum of Vq corresponds to the uniform 
configuration A(x) — Aq at a = without a split-off level. The function Vo(a) monotonically increases from Vb(0) = 
to Vo(oo) = 2W S , ( W s — 2/7tA is the total soliton energy) with the asymptotic behavior. 

V (a) w — A a 3 . a -> 0, V (a) « 2W S - 2E (a), a oo. (34) 

For the polaron v = 1 (either one electron at Eq or one hole at —Eq) the equilibrium state corresponding to the 
minimum of V is achieved at a = ao where 

sinha = 1 , mmVx= V(a Q ) = W p = 2 3 / 2 A /vr, E = A /v^. 

The limiting values of Vi(a) are Vi(0) = A , Vi(oo) = 2W S . 

The case v = 2 corresponds to either an exciton (one electron and one whole at ±E that is each of these levels 
is singly occupied) or to electron or hole bipolarons (all states ±£-are either filled or empty). The excitonic state 
plays a principal role in the subgap optical absorption problemEJEl The equilibrium state is achieved at a = oo 
where V — 2W S and Eq = 0. The function V2(a) decreases monotonically from V2(0) = 2Ao to V<j(oo) = 2W S . The 
dependencies V u {a) arc shown in Fig. 3. 
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V 




FIG. 3. Plots of Vv(a) for v = 0, 1, 2. 



The above a— dependent family of solutions is the only one which provides exactly one pair of discrete levels in 
the gap. Its perturbations originate only shallow levels located near the edges ±Ao which are not important for the 
described processes. Therefore we will choose the configuration A(x,t) = A s (x,a(t)) of ( |32] ) taken with the time 
dependent parameter a. The kinetic term in the Lagrangian becomes f(a)d 2 with 



/«,<»)) = (1= 



(35) 



The function /(a) increases from /(0) = reaching the maximum value / ~ 1.55Ao/(7tAojq) at a w 0.78, and then 
decreases with asymptotics 



/(a) 



2 Co Aoa 

3 TT UJq 



0: 



/(ao) 



1.53 A 



/(oo) 



2 A 

3 7TWq 



with the same coefficient Co as in the previous section. 
The translational mass becomes 



(36) 



M(a) 



tlx 



dA 

dx 



8Aj_ 

9 2 ^o 



tan a a cosh a — sinh a 



cosh 3 ( 



The function M(a) monotonically increases with asymptotics 

0; M (ao) 



M(a)«C M ^|a 5 , C M 



32 
15' 



A 3 

0.49^ 



M(oo) 



16 Af 

37T 



(37) 



(38) 



Following the prescriptions of Sec. II, we can describe now the PG region W p < SI < Aq. In the limiting cases the 
results are qualitatively similar to the ones for the shallow polaron. 
Near the free edge, SI w Aq we find 



i (si) ~ (A - «r 1/4 



exp 



32VC0" (A - si) 3 / 2 



9tt 



(39) 



with the same coefficient Co as in (p9|). 
Near the polaronic energy SI « W p we find 



/(il) = const' 



1 



exp 



-Ci 



Ao 

w 



exp 



m-Wp) C3A0 

C 2 — log- 



w 



(fi - W p ) 



(40) 



Here the numerical coefficients C\ = 0.34, C2 = 5.2, C3 =0.1. 
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D. ARPES intensities. 



For case of ARPES the potential V\ acquires an additional term P 2 /2M, therefore the closed trajectory is allowed 
in the region (fl, P) where (See Fig. 4) 

Vi(a ro ) = Vi{a m ) + P 2 /2M(a rn ) - fi < , dVi/da| a=am = (41) 

Here the value a m is the point of a minimum of the potential Vi(a), which exists for any value of the momentum P. 
In the limit P — we have a m = oq. For other f2, P there are no solutions with a finite action so that the transition 
rate is zero. In contrast to the PES case, the closed trajectory exists for any frequency f2 > W p in the region of the 
momentum P satisfying the condition (p4|). 




FIG. 4. Plots of V (a) and Vi + P 2 /2M(a) - fi. 



Consider most important cases where analytical solution can be found. 

A. In a vicinity of the absolute edge fi w W p the absorption appears at £1 > V\(a m ) + P 2 /2M(a m ) = W p . Here, 
near the polaron edge, the translational mass is almost constant staying near M(a m ) hence the momentum appears 
simply through the shift Q =>• O — P 2 /2M(a m ). 



i(p, n) oc 



n-w„ 



:exp[-5 ]. 



(42) 



with 



So 



Ci Ap 
s/ttX w o 



G 



C 3 A C 



ttXloq (fl — W p ) 



(43) 



The result for the absorption looks like the Eq. (| 
functions of P: C\, C-x, C3. They are defined as 



for the PES with W p =>■ W P (P) but the coefficients Cj become 



d = 4VttA 



A 



da^f(a)V (a)+ / daJ f(a)(Vi - W p ) 



C 2 = VttXloq ; , C 3 = -(df - a T ) — -r . 

V»(a m ) 2 Ao 

and at P — > they coincide with the numerical ones Ci, C2, C3 defined above in Sec. IIIC. 
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B. A vicinity of the free edge fl w Ao is very particular for the ARPES. It is dominated by shallow fluctuations 
determined by small a <C 1 where the kinetic energy diverges as P 2 /M ~ a~ 5 according to (|38l). The form-factor 
| | 2 acquires an explicit form 

I^a^" 2 ^ < 44 > 

To simplify appearance of cumbersome relations we shall use below dimensionless units Ao = 1 for the momentum P 
and energies which also will be counted with respect to Ao. In the limit et t , a/ <§C 1 equations ([[70, ( |lS| ) are reduced 
to 

-^-T + ^r 83 ' (45) 

2 GmOIt 



3tt 



4 



2 



0, 



(46) 



where e = (fi — A )/A and the term comes from the expansion of the local level energy E . A simple analysis 

of (f45|), (p6|) shows that af — a? = o(a^), and first term gives the main contribution to the action (pOj), so that 



S = 



9tt ojq 



(47) 



1/3 

The general constraint is So S> 1 hence ar ^> w ■ ^ assures us also in the condition M ^> Ao (the mass M must 
stay above the free band edge mass) which requires for a weaker inequality a-r^^- 

We have different regions where one of terms of Eq. ( ^B| ) is small in comparison with others. 

Bl. Consider at first the region lying deeply enough within the PG: 

{uJoP) i/7 < -e < 1 

Here the characteristic values of ax are not too small so that one can neglect the kinetic energy ~ P 2 ^o l a T m compare 



to the energy — a^/2 of the localized level: P uj^ja^ <c ai. Then ot 



/u>o and d£l/ dT ~ ot- The 

momentum dependence comes only from the form-factor ([44j) but it can be still appreciable: P/ax is limited by a big 

quantity (P/cit <C e 5 / 4 /w ^ w o tnat iS allowed to be large. 
The total intensity is 



1 3/4 



exp 



32\/G^ 

97TW0 



|3/2 



ttP 



2e IVa 



(48) 



In the limit — e <C \/Pujo the contribution from | ^> v | 2 (the second term in the exponent) dominates, while in the 
opposite case the main contribution comes from the term So- The line of the maximum intensity — e ~ V^qP will 
show up in the ARPES plots as a quasi spectrum which intensity decreases exponentially with growing momentum: 
I(P) ~ P- 3 / 8 exphconstP 3 /^" 174 ]. 

B2. Consider a marginal frequency lying very close to the free edge e = for both signs of e: 

uj q < |e| < {uJ a P) i/7 < 1 for P > uj 2/5 . 
In this case the instanton kinetic energy and the binding of the electron almost compensate each other: P 2 /2M ss 

1 1 7 

a\j2 3> |e| which gives us qt = {[ujqP) 2 / C m) and we obtain 



7(e,P)oc 



p3/7 



exp 



71 2/7 r 



(49) 



The form-factor ^> 2 p gives the main contribution (the first term) to the exponent of 
B3. Consider now a new region of positive frequencies e > 0: 



» (w P) 4/7 but e « ^ /3 P 2 . 
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The first inequality ensures that the excess frequency is absorbed mainly by the kinetic energy with a negligible 
(afn/2 <C e) contribution from the electronic state which is still localized. The second inequality supports the adiabatic 
condition So ^> 1. We find 

aT ~ (2C M e) 1 /5' dtt^T^' 

The ARPES intensity becomes 

pl/S 

J ( £ ' P ) a ^io ex P 

The form-factor contribution, the first term in the exponent, always dominates but the second term is also big thus 
contributing to the dependence J(e, P). The total expression in the exponent of (|5^) looks non monotonous similar 
to the case Bl but now the minimal is not physical: its position e ~ (ujqP) 3 ^ 4 would fall too low, into the region B2. 
Contours of constant intensity are close to the descending line e ~ w 2 /P 3 . The contribution ( |50| ) from the fast moving 
instantons can be observable only below the intensive free electronic absorption at e w P 2 /2 that is at P ^ lu^ 5 (a 
similar constraint already appeared for the case B2). 

For the ARPES intensities we must, in principle, minimize the total effective action with the prefactor term | iffp | 2 , 
since in some cases this term gives a main contribution to the exponent. Then the action is changed as 

S^S-log | * P | 2 

The extremal solution, as before, satisfies the Eq. (|p]). The turning point a/ is not changed and is defined by the 
Eq. fllSj), but the equation defining the point ax becomes more complicated. Instead of the Eq. (|l7]) we obtain 

(51) 

a—aT 

This equation leads to the discontinuity of the potential V at the point ar and, due to the equation of motion ( |1S|) , 
leads also to the discontinuity of the "velocity" a. The expression for the ARPES intensity has the same form ( |19| ) 
with the action (p0|). The only difference is the shift of the point ot due to Eq. (pTf). 

For a considered in B1-B3 regions the point ar is defined by ( fl5| ) with the additional term ~ Puio/ar- The simple 
analysis shows that this term is small in comparison with the kinetic part P 2 /2M or with the term a 2 /2, therefore it 
can be neglected. Thus, the main contribution to the action is not changed, and all above results are retained. 



ir(2C 



M 
2/5 



,1/5 



p3/5 £ l/5 



9tt(2Cm) 3/5 e 3 / 5 



(50) 



~ l gl0gj_gi 

V {a T ) - Vi(a T ) = ^ — 




E. Optical absorption: effects of confinement. 

The same method is applicable to the subgap optical absorption problem which was basically studied already i 1100. 
We shall see that more details can be easily obtained. In this case we must use the potential Va, defined in (|3^), 
instead of V\. The term V2 describes the first excited state when the optical photon, assisted by a lattice quantum 
fluctuation, creates an e — ft, pair with levels ±E(t) spanning the whole interval \E\ < Aq. The transition rate is given 
by the modified Eq. (^). For the vicinity of the free electronic edge f2 ~ 2Ao we easily find, in analogy with (|3S|), 
that 



oc exp 



Q) 



3/2 



9tt 



(52) 



This law was noticed in@, while with a differerLt-jCocfficient in the exponent. Actually it was obtained already inBby 
another method which was later reproduced intH. Near the absolute absorption edge Q w 2W S — 4Aq/tt we obtain 



oc exp 



(53) 



where /(oo) is given by Eq. (p6|). Near this edge the absorption is due to tunneling to the final state formed by the 
two diverging kinks with the distance parameter a — > 00. (This regime was studied for the first time inE3 but the 
above analytical asymptotics was not found.) The law (p3h is different from (E^) which was typical for all polaronic 
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thresholds. The particularity is related to a very fast, ~ expj— a], decrease of interactions for diverging solitons which 
determine the threshold. In terms of the effective particles the difference is that near the polaronic threshold the 
turning point approaches the potential extremum while for the solitonic threshold it climbs asymptotically to the 
plateau (the dashed line at the Fig. 3). 

While the law @ is very robust, e.g. the same as for shallow polarons, the solitonic threshold is extremely sensitive 
to perturbations. The most drastic effect comes from the confinement energy 8V U = Fa with the constant confinement 
force F originated by lifting of the ground state degeneracy. The confinement is always produced by the interchain 
couplingtj. But even for a single chain it appears ia. .cases of a build-in alternation of unit cells which interferes with 
the spontaneous dimerization like in cis — (CH)xc3'£J. In principle, the solitonic terms for this problem have been 
found exactly in the frame of the microscopic model of the "combined Peierls state"EaO and they can be used in our 
calculations for an arbitrary strength of the force F = 7A /£o- But here we shall consider only a weak confinement 
7 <C 1 which preserves the local structure of solitons but prevents their divergence at large a. This effect can be taken 
into account by adding the term 7A a to the potentials Vb, Vi taken in zero order in 7 -C 1. (For the PES this effect 
is not relevant since the minimum of the polaronic potential V\ is achieved at finite value do where confinement is not 
important yet.) As a result the minimum of the potential V2 is shifted to the finite point a 7 log(16/7T7)/2, where 
V 2 "(a 7 ) = 27. The optical absorption edge is shifted up to W 1 = 4Ao/7r + 7A0 log(16/7T7)/2 and the transition rate 
1(0,) is given by Eq. @ if substitute W p W 7 , VI -> V2, a Q -> a 7 , a x -> cosh"ii^/2). 

The problem of solitonic absorption, without the confinement, was considered inEj'tj by means of a different quasi- 
classical approaches. Our results supply the absent analytical asymptotics in this case and provide a new description 
for the general case of confinement. Further complications are coming from the long range Coulomb interactions 
which must be studied before comparison with existing experiments. 



F. Quantum fluctuations as an instantaneous disorder. 



In a one dimensional system the X-ray or optical absorption near the band edge can be viewed as the one in a 
system with a quenched disorder which is emulated by instantaneous quantum fluctuations. This conclusion was 
achievedlij by analysis of the perturbation series for the polarization operator which is not the most efficient way to 
do. Later applications of this idea did not add any prove£3. Here we shall derive this result within the approach of 
functional integrals which will provide a model independent limit for over approximate analysis of more complicated 
regimes and systems. 

Consider a particle interacting with a media of harmonic oscillators. This model embraces the shallow polaron case 
of Sec.B and the limit of Vt « Aq of Sec.D. The probability is given by the action 



S{Q,^-T) 



dPx 



+ / d°x 



dt 



dt 



dQ 

dt 



1 <9* 

2m dx 



(54) 
(55) 



We can integrate out the field Q at all x and t to arrive at the action which is defined only at the interval (0, T), c.f.t 



S{^;T) 



d u x 




dx 



1 2 
p wo 



dt\ / dt2p{x,ti)p(x,t 2 )ex.p[-uj \ti 



where p = | \T/ 1 2 . The result will prove that the characteristic energy scale Sd is much larger than ujq which takes 
place only at D < 2. Then the whole time interval is short Two ~ wq/5d *C 1 and we can neglect the retarda- 
tion: exp[— uJo\ti — ^2 1] — ► 1- Now the last term in S{^;T) can be decoupled back by the Hubbard-Stratonovich 
transformation via the time independent auxiliary field £ to give us 



won-/< D .[jf*(s : |£r) + «^ 



dt 1 p(x,t) + ^—( 2 (x) 
g z uj 



Finally we can integrate over ^ to arrive at the transition probability 

1 



dTexp 



g 2 uj 



?(x)+T(E[ax)]-Sl) 
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After rotation to the real time it becomes 



D[((x)S(E[((x)} - ft) exp 



where is the eigenfunction in the random field £ 



ff 2 ^o 



C 2 0r) 



Clearly the dispersion <5o = S-y/wo/2 of the field £ is just the mean square-pf quantum fluctuations of the phonon 
field gQ(x,t). The well known results for the DOS of disordered systemsErLil confirm our direct calculations for the 
time dependent processes in the pseudogap region and extend them to the nonexponential part of the spectrum of 
CI > 0. For D = (m = oo) and D — 1 we arrive correspondingly at the widths So ~ and Si ~ (<? 2 cc>o) 2 / 3 to 

find that in both cases the condition 5u ujq is satisfied at low enough u>q. (This condition is not satisfied at higher 
dimensions hence at D > 2 the quantum fluctuations are not reduced to the random static field but rather become 
resolved phonon assistant processes.E3) An additional peculiarity of the Peierls model is that, in microscopic units of 
A and £o the coupling constant g = 1, then the band edge smearing is given by Si = (Aqw 2 ) 1 ^ 3 . 



IV. DISCUSSION AND CONCLUSIONS 

Electronic properties of quasi one - dimensional conductors in the Peierls state are quite peculiar in several respects. 
In general it is due to a strong interaction of CDW deformations with normal electrons which leads to their fast self 
trapping. The stationary excited states of the model are solitons (kinks) and polarons with energies E s ,E p < Ao. 
The processes related to these nonlinear excitations determine the true gaps 2E s ,2E p placed within the pseudogap 
2Ao. The single particle gap, as measured in absorption or tunneling, is opened by polarons which should exist also 
in ID systems without symmetry breaking, like the majority of conjugated polymerstjO. The effect is very common 
because in ID semiconductors the selftrapping of free electrons takes place for any type of e-ph interactions^ 2 ] while 
in higher dimensions the long range interactions are required. Here the minimal role of the Peierls effect, or at least 
its partial contribution, is to ensure a presence of the strong e-ph coupling. The optical threshold exists at 2E p 
in general. But for systems with degenerate ground states like the ■ pplyacethylene there is also a lower gap at 2E S 
because the lowest excitations are now topological solitons (kinks )EjSE3. A weak interchain coupling^ preserves the 
polaronic effects, while creating a shallow barrier with respect to the selftrapping, thus allowing for metastable free 
electronic states and corre sponding transitions. The 31? effects on solitonic transitions are more drastic because of 
the confinement effectE§E3. 

Presented results recover the stretch exponential dependencies near the PG edge (for both PES and optics) 



exp 



ci 



(56) 



with different powers !/£>>!. Near the absolute threshold at SCI = Cl — W p <C W p we find the same (for the PES) law 



exp 



UJ 



UJQ oi i 



(57) 



which differ only by numerical coefficients among different models. The first constant term in the exponent describes 
the adiabatic reduction for the probability of creation of the stationary polaron at W p . The second term gives the law 

for approaching this threshold. At first sight, this quasi linear dependence of the exponent ~ £p In is weaker than 



the one ~ In coming from the pre-exponential factors ~ (SCI) 11 but it is not true: by definition we are working 

at SCI which can be small only in compare to W p but must be large in compare to uo. Hence ^ In -^y 3> In ^ and 
the exponential term determines the profile near the absolute threshold. The same law holds for the two particle 
process (the internal optical absorption) if the ground state is nondegenerate. But for a system with a spontaneous 
symmetry breaking like the Peierls model for the polyacethylene the threshold dependence of the exponent changes 
from ~ SClhiSCl to a stronger one ~ \/Cl — 2W S with the threshold 2W S being the energy of a solitonic pair. 

The momentum dependence of the intensity I(Cl, P) as recovered by the ARPES shows a rich variety of regimes. 
Only near the absolute threshold the law (p>7]) is simply generalized by adding the polaron kinetic energy: W p => 

Eg demonstrates several nontrivial regimes. The 



W p + P 2 /2M p . But the region near the free electronic threshold Cl 
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differences come from two effects achievabie for shallow states. One is the momentum dependence of the form factor - 
the momentum distribution of the split off state. Another one is the inertial part of the action coming from a drastic 
dependence of the kinetic mass on the localization parameter which evolves in time along the instanton process. One 
of unexpected results is that at large P the law (|56| ) changes to a non monotonous function of VL which would show 
itself as a quasi spectrum (a line of maxima of P)) within the pseudogap. Another result is an appearance of 
optimal localized fluctuations at an elevated kinetic energy so that they show up above the PG region at f2 > 0. 
These observations may be more general than the studied one dimensional adiabatic models. 

In conclusion, we have calculated intensities of subgap photo electronic transitions by means of a functional in- 
tegration over the lattice oscillation modes. By this virtue we have studied the problem of photoemission from the 
pseudogap region for typical one-dimensional models. We have found general expressions for the transition rate and 
studied the asymptotic behaviors near absorption edges, both below the free electronic edge and approaching the lower 
true one. It was shown that the main contributions to the transition rates comes from instanton configurations of 
the phonon field. The Peierls model was investigated in particular details in relation with its numerous applications. 
Peculiarities of the ARPES regime are coming from unexpected effects of the instanton motion. 

For the intragap optical absorption problem we have calculated the asymptotic behavior of the absorption coefficient 
near the threshold for creation of pairs of solitons and studied effects of confinement. 
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